A homoclinic tangle on the edge of shear turbulence 
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Experiments and simulations lend mounting evidence for the edge state hypothesis on subcritical 
transition to turbulence, which asserts that simple states of fluid motion mediate between laminar 
and turbulent shear flow as their stable manifolds separate the two in state space. In this Letter we 
describe a flow homoclinic to a time-periodic edge state. Its existence explains turbulent bursting 
through the classical Smale-Birkhoff theorem. During a burst, vortical structures and the associated 
energy dissipation are highly localized near the wall, in contrast to the familiar regeneration cycle. 

PACS numbers: 47.27.Cn, 47.27.ed, 47.52.+j 
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Introduction In recent years, the open problem of 
subcritical transition to turbulence in shear flows has 
seen a surge of interest, sparked by a rapidly increas- 
ing ability to perform numerical simulations as well as 
by a string of novel applications of dynamical systems 
theory. The dynamical systems approach advocates the 
idea, that transitional shear flow is regulated by special 
solutions with a relatively simple spatial structure. They 
may be travelling waves, time-periodic solutions or even 
solutions chaotic in time. Such states live in between a 
stable laminar state and stable, or met a-st able, turbu- 
lence. In simulations, these special states can be iden- 
tified by a shooting algorithm in which initial data are 
iteratively refined to yield a flow which neither laminar- 
izes nor becomes turbulent, but instead lingers on the 
"laminar-turbulent boundary" [1]. States on this bound- 
ary are necessarily unstable and consequently they can 
only be observed as transient effects in experiments |2i]. 
Their relevance to the transition process lies mainly in 
their stable and unstable manifolds in state space. These 
determine how the fluid behaves as it transitions from 
near-laminar to turbulent states and vice versa. 
A particularly interesting situation arises when the 
laminar-turbulent boundary is formed - at least locally - 
by the stable manifold of a traveling wave or periodic so- 
lution, which is then called an "edge state" . Edge states 
have now been computed for flow in channels as well as 
pipes, with various numerical schemes and spatial dis- 
cretizations, and there is growing consensus that they 
are a robust feature of subcritical shear flow [3] . 
Logically, the next step in this analysis would be to study 
the stable manifolds of the edge states. Knowledge of 
the geometry of these manifolds would open the door to 
the application of control techniques, aiming at a forced 
laminarization of the flow. Indeed, some results to this 
effect have been obtained using linearization about an 
edge state Little, if anything, is know about the 



global, nonlinear structure of these separating manifolds, 
greatly reducing the usefulness and predictive power of 
the edge state hypothesis. The direct study of a separat- 
ing manifold is hard if not impossible, owing to its high 
(formally infinite) dimensionality. In the current Letter, 
we commence by studying the two-dimensional unstable 
manifold of an edge state in plane Couette flow, using a 
novel computational algorithm. We find that it contains 
an orbit which returns to the edge state along its stable 
manifold. Through the classical Smale-Birkhof theorem 
|5j, the presence of this homoclinic orbit implies the ex- 
istence of an intricate tangle of the stable and unstable 
manifolds and chaotic dynamics which manifests itself as 
irregular turbulent bursting. 

Transitional plane Couette flow We consider plane 
Couette flow at a Reynolds number of Re = 400 in the 
minimal flow unit of dimensions L x W x H ^ where Re 
is based on half the velocity difference between the two 
walls, /7/2, and half the wall separation, H/2. The 
streamwise and spanwise periods are {L/H^W/H) = 
(2.76,1.88) [6]. We used resolutions of 16 x 33 x 16 as 
well as 32 X 33 X 32 grid points in the streamwise (x), 
wall-normal {y) and spanwise (z) directions, respectively, 
and checked that the behaviour is qualitatively the same. 
In this computational domain, no Nagata steady so- 
lution exists [8]. Instead, the edge state in this flow is a 
time-periodic variation of the laminar flow, which shows 
weak, meandering streamwise streaks 0, Q. This gen- 
tle Unstable Periodic Orbit (UPO) has a single unsta- 
ble Floquet multiplier and thus its unstable manifold has 
dimension two and its stable manifold has codimension 
one in state space. Both the UPO and its unstable ma- 
nifold are contained in a subspace invariant under the 
spatial symmetries given by reflection in the mid plane, 
followed by a streamwise shift over L/2, and reflection in 
the streamwise and spanwise direction, follow by a span- 
wise shift over W/2. Consequently, we can impose these 
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symmetries on the solutions to reduce the number of de- 
grees of freedom, without placing artificial restrictions on 
the fluid motion. At the higher resolution, the number 
of degrees of freedom in the simulations is about 11, 000. 
Homoclinic orbit computation The unstable mani- 
fold of the UPO can be computed using multiple-shooting 
orbit continuation [lOj. Essentially, this algorithm pro- 
duces a sequence of orbits, contained in the manifold, 
by arclength continuation of a boundary value problem 
in time. Because of the extremely sensitive dependence 
on initial conditions in turbulent Navier-Stokes flow, we 
need to compose the orbits of multiple segments, each 
not much longer than the decorrelation time, i.e. the 
time scale of exponential divergence of turbulent states 
of fluid motion which are initially close. The results pre- 
sented below use up to six shooting intervals, the inte- 
gration time on each interval staying below two times the 
period T of the UPO and below five times the decorrela- 
tion time. 

As we compute a sequence of orbits, it may happen that 
it converges to a homoclinic orbit, which separates the 
unstable manifold into two components. Such an orbit is 
shown in Fig. [Utop). In this projection on energy input 
and dissipation rate, the UPO is the tiny loop labelled 
L. In the background we have plotted the Probability 
Density Function (PDF) of transient turbulence. In the 
transition to turbulence the homoclinic orbit overshoots 
the maximum of the PDF, then passes close to it on the 
way back to the UPO. Clearly, the shape of the PDF 
agrees well with that of the homoclinic. 
We generically expect the homoclinic orbit to approach 
the UPO along the least stable subspace of the latter, 
also called the leading stable subspace. This is visibly 
the case in FigfT^bottom), which shows a close-up of the 
UPO, along with the local unstable manifold in red and 
the local leading stable subspace in blue, now with a 
third axis showing the mean square streamwise vorticity. 
How well the depicted orbit approximates a homoclinic 
connection can be quantified by measuring two distances: 
one from the end point of the computed orbit to the lead- 
ing stable subspace at a point on the UPO and one to 
this point along the leading stable subspace. We found 
these distances to be of order 10~^ and 10~^ in the en- 
ergy measure, respectively. To test the robustness of this 
result, we have recomputed the homoclinic with a smaller 
integration time step and a varying number of shooting 
intervals. 

The flow structure of bursting Figure [2] shows the 
time-evolution of flow structures along the homoclinic or- 
bit at the six phases indicated in Fig.[Tl In the early stage 
of the evolution (phase a) the spanwise standing-wave 
motion of the streak is enhanced, so that the stream- 
wise dependence, i.e. the three-dimensionality, of the 
streak gradually becomes significant. At the same time 
the streak grows in the wall-normal direction. Such be- 
haviour is a direct consequence of the linear instability of 
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FIG. 1. Top: projection of the homoclinic orbit onto energy 
input and dissipation rate, normalized by their value in lam- 
inar flow, e and e. The piece of orbit leaving the edge state, 
L, is shown in red and the one approaching it in blue. In the 
background, the PDF of transient turbulence is shown in gray 
scale. The labels a-f correspond to the snap shots in Fig|2] 
The inset shows the input (red) and dissipation (green) rate of 
energy as a function of time, normalized with T, along part 
of the homoclinic orbit. Bottom: close-up of the UPO, L, 
with the homoclinic. The red and the blue arrow denote the 
projected local unstable manifold and the projected leading 
stable subspace, respectively. On the axes are the normal- 
ized energy input and dissipation rate and the mean square 
streamwise vorticity normalized by the mean square total vor- 
ticity, where the mean is taken over the flow unit. 



the UPO. The eigenstructure for the instability is char- 
acterized by disturbances of the streamwise velocity and 
vorticity which are highly localized on the crest and the 
valley of the streak. 

As time progresses, the simultaneous spanwise oscillation 
and wall-normal growth of the streak exceed a critical 




FIG. 2. Visualizations of flow structures in one periodic box L x W x H in six phases on the homoclinic orbit, labeled as in 
Fig.[T] In each phase the structures are viewed in two distinct directions. Gray corrugated isosurfaces of the null streamwise (x) 
velocity represent streamwise streaks. Red and blue objects are isosurfaces at 0.2(11/ H)"^ for the second invariant of a velocity 
gradient tensor, and denote the vortex tubes of the positive (clockwise) and negative (counter-clockwise) streamwise-vorticity 
component. Green isosurfaces show the local energy dissipation rate at 20 times the value e in laminar flow. Cross- stream 
{y,z) velocity is shown in the mid plane x — L/2. By the flow symmetry it is related to that at x = by a reflection in the 
spanwise direction. 



level for bursting. There appear thin layers of extremely 
high vorticity between the wall and the crest and the val- 
ley of the streak, in which intense quasi-streamwise vor- 
tex tubes are generated to induce high shear and thereby 
high energy dissipation (phase b). The quasi-streamwise 
vortices of opposite signs of the streamwise vorticity align 
following the oscillation of the streak in the spanwise di- 
rection. Then the streak is rapidly deformed, while the 



second pair of quasi-streamwise vortices appears on the 
crest and valley, and align with the first pair to form an 
array which moves fast on the crest or valley in the span- 
wise direction (phase c). During this streak deformation 
and vortex generation intense energy dissipation is ob- 
served in the thin layers between the wall and the crest 
or the valley of the streak. Streamwise vortices similar 
to those found in the regeneration cycle [6] and the cor- 
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FIG. 3. Temporal variation of energy-dissipation (solid 
curve) and volume (dotted curve) fractions in the region 
bounded by the green surfaces in Fig. [21 where the local dissi- 
pation rate is greater than 20 times the value in laminar flow. 
Time t is normalized with the period of the gentle UPO, and 
the six phases on the orbit are shown by dashed vertical lines. 



responding strong UPO [9] also appear on the flanks of 
the streak. However, the streamwise vortices on the crest 
and valley of the streak are not observed in the near- wall 
regeneration cycle, and therefore we can say that these 
vortices are typical flow structure of bursting with signif- 
icant energy dissipation. Actually, an inspection of the 
turbulent state of plane Couette flow has shown that in 
bursting events of intense dissipation there appear similar 
vortex arrays on the crest and valley of the highly grown 
streak, which are associated with strong local dissipation. 
As shown in Fig. [H the intense local dissipation (repre- 
sented by green isosurfaces) observed around the crest 
and valley of the streak in Fig. [2] (phases c and d) con- 
tributes to the rapid increase of total energy dissipation 
(see Fig. [T]). The identified high-dissipation regions oc- 
cupy less than 10% of the spatial domain but account for 
over 40% of the energy dissipation, implying that strong 
energy dissipation in the bursting event can be attributed 
to flow structures localized in the region around the crest 
and the valley of the streak shown by the green surfaces 
in Fig. [21 As time goes on, the highly deformed streaks 
are broken down (phase d) , and then the streak and vor- 
tices decay rapidly (phase e) as the flow finally returns 
to a quiescent state close to the gentle UPO (phase f). 
Transient approaches of a turbulence state to relatively 
quiescent states of low energy input and dissipation are 
followed by bursting events with intense energy dissi- 
pation, as reported by Kawahara & Kida [9], who sug- 
gested the interpretation of bursting in terms of a two- 
way heteroclinic connection between the gentle UPO and 
a UPO embedded in turbulent flow. The present more 
rigorous analysis has presented another interpretation. 



namely that the whole process of the bursting, including 
a transient approach to a quiescent state and a subse- 
quent highly dissipative event, is the manifestation of a 
homoclinic orbit arising from the gentle UPO alone. 
Conclusion The existence of an orbit homoclinic to 
this edge state implies that the geometry of the laminar- 
turbulent boundary is rather complex and we can expect 
to find a manageable approximation to it only locally. 
At the same time, it implies the existence of infinitely 
many UPOs which correspond to flows with arbitrarily 
many, arbitrarily long, near laminar izat ion events. It is 
natural then to think of turbulent shear flow as governed 
by a large chaotic attractor which comprises both pe- 
riodic orbits in the turbulent regime, which reproduce 
the regeneration cycle and periodic orbits which re- 
produce near-laminarization and bursting events. Phys- 
ically, these bursting events are very different from the 
regeneration cycle. The streamwise vortices appear on 
the crest and valley of the streaks rather than on their 
flanks, and the larger part of the energy dissipation takes 
place in the cycle around the vortices in the near-wall, 
high strain region. 

The homoclinic solution presented here adds a new el- 
ement to the elucidation of turbulence by means of dy- 
namical systems theory, namely that of temporal local- 
ization. Since it has been computed in a minimal flow 
unit, in cannot capture spatial intermittency. Recently, 
several equilibrium and travelling wave solutions exhibit- 
ing spatial localization have been found [llj . A clear goal 
for the near future is to find periodic or connecting orbits 
which combine the two forms of localization. 
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